Differential Allele-specific Expression Research Articles

Single-cell allele-specific expression analysis reveals dynamic and cell-type-specific regulatory effects

Differential allele-specific expression (ASE) is a powerful tool to study context-specific cis-regulation of gene expression. Such effects can reflect the interaction between genetic or epigenetic factors and a measured context or condition. Single-cell RNA sequencing (scRNA-seq) allows the measurement of ASE at individual-cell resolution, but there is a lack of statistical methods to analyze such data. We present Differential Allelic Expression using Single-Cell data (DAESC), a powerful method for differential ASE analysis using scRNA-seq from multiple individuals, with statistical behavior confirmed through simulation. DAESC accounts for non-independence between cells from the same individual and incorporates implicit haplotype phasing. Application to data from 105 induced pluripotent stem cell (iPSC) lines identifies 657 genes dynamically regulated during endoderm differentiation, with enrichment for changes in chromatin state. Application to a type-2 diabetes dataset identifies several differentially regulated genes between patients and controls in pancreatic endocrine cells. DAESC is a powerful method for single-cell ASE analysis and can uncover novel insights on gene regulation.

Differential Allele-Specific Expression Revealed Functional Variants and Candidate Genes Related to Meat Quality Traits in B. indicus Muscle

Traditional transcriptomics approaches have been used to identify candidate genes affecting economically important livestock traits. Regulatory variants affecting these traits, however, remain under covered. Genomic regions showing allele-specific expression (ASE) are under the effect of cis-regulatory variants, being useful for improving the accuracy of genomic selection models. Taking advantage of the better of these two methods, we investigated single nucleotide polymorphisms (SNPs) in regions showing differential ASE (DASE SNPs) between contrasting groups for beef quality traits. For these analyses, we used RNA sequencing data, imputed genotypes and genomic estimated breeding values of muscle-related traits from 190 Nelore (Bos indicus) steers. We selected 40 contrasting unrelated samples for the analysis (N = 20 animals per contrasting group) and used a beta-binomial model to identify ASE SNPs in only one group (i.e., DASE SNPs). We found 1479 DASE SNPs (FDR ≤ 0.05) associated with 55 beef-quality traits. Most DASE genes were involved with tenderness and muscle homeostasis, presenting a co-expression module enriched for the protein ubiquitination process. The results overlapped with epigenetics and phenotype-associated data, suggesting that DASE SNPs are potentially linked to cis-regulatory variants affecting simultaneously the transcription and phenotype through chromatin state modulation.

Abstract 5725: HOXA9 demonstrates allelic imbalance in KMT2A-rearranged infant acute lymphoblastic leukemia

Abstract Background: KMT2A-rearrangement is a characteristic finding in approximately 70% of infant acute lymphoblastic leukemia (ALL) cases. We previously utilized a computational screen to find that KMT2A-r infant ALL demonstrated somatic allele specific expression for 431 genes when compared to non-KMT2A-r samples (Abstract Number 3649, AACR 2019). The gene that exhibited the most differential allele specific expression (ASE) pattern was HOXA9. Here, we follow up that analysis with manual examination of additional infant ALL samples. Methods: We performed whole genome sequencing (WGS) and RNA sequencing (RNAseq) on peripheral blood or bone marrow samples from 48 infants with ALL (33 KMT2A-r cases and 15 non-KMT2A-r cases) treated on Children’s Oncology Group protocol AALL0631, at diagnosis (DX), at remission (MD), and at relapse (RL), if applicable. WGS data at MD was phased using a 1000 Genomes reference panel. HOXA9 expression was quantified in transcripts per million (TPM). Biallelic expression was manually calculated on the phased genome for samples with at least 1 germline single nucleotide polymorphism (SNP) within HOXA9 and at least 5 aligned RNAseq reads over the HOXA9 transcript. ASE was defined as minor allele fraction (MAF; allele fraction for the less prevalent allele) of ≤ 0.35. Whole genome bisulfite sequencing data and sequencing-derived copy number data were also analyzed to investigate possible explanations for HOXA9 ASE. Results: We identified diagnostic samples from 21 patients, 9 KMT2A-AFF1, 9 KMT2A-MLLT1, and 3 non-KMT2A-r, that passed manual filtering. Data included 5 patients with paired DX and RL samples. Of the 15 KMT2A-r samples unable to be analyzed for ASE, 13 lacked a SNP for phasing and 2 had < 5 reads. For the 12 non-KMT2A-r samples unanalyzable for ASE, 6 lacked a SNP and 6 had < 5 reads. Average HOXA9 expression was 15.17 TPM (range 0.01 to 44.34) in KMT2A-r samples and 2.91 TPM (range 0 to 29.29) in non-KMT2A-r samples. ASE of HOXA9 (MAF ≤ 0.35) was observed in 14 samples, 72% of KMT2A-r (13 of 18 samples) and 33% of non-KMT2A-r (1 of 3 samples). RL samples were found to have HOXA9 ASE in all cases with ASE of HOXA9 in the DX sample. Of the 18 KMT2A-r samples, 6 of 9 with KMT2A-AFF1 and 7 of 9 with KMT2A-MLLT1 demonstrated HOXA9 ASE. Frequency of HOXA9 ASE at diagnosis was similar between patients known to later relapse versus those known to remain in remission. Allele specific methylation, mono-allelic deletion, or copy-neutral loss of heterozygosity (cn-LOH) for HOXA9 were not observed in KMT2A-r nor non-KMT2A-r samples. Conclusion: Our expanded cohort supports the ASE of HOXA9 in KMT2A-r infant ALL. This phenomenon did not appear to segregate by KMT2A-r partner gene, nor could it be explained by allele specific methylation or copy number variation (e.g., deletion of one HOXA9 allele or cn-LOH). Overall, further investigation into a causal mechanism is needed. Citation Format: Sarah E. Mc Dermott, Byunggil Yoo, Warren A. Cheung, Emily Farrow, Bing Ge, Margaret Gibson, Patrick A. Brown, Erin Guest, Tomi Pastinen, Midhat S. Farooqi. HOXA9 demonstrates allelic imbalance in KMT2A-rearranged infant acute lymphoblastic leukemia [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2022; 2022 Apr 8-13. Philadelphia (PA): AACR; Cancer Res 2022;82(12_Suppl):Abstract nr 5725.

Improved haplotype inference by exploiting long-range linking and allelic imbalance in RNA-seq datasets

Haplotype reconstruction of distant genetic variants remains an unsolved problem due to the short-read length of common sequencing data. Here, we introduce HapTree-X, a probabilistic framework that utilizes latent long-range information to reconstruct unspecified haplotypes in diploid and polyploid organisms. It introduces the observation that differential allele-specific expression can link genetic variants from the same physical chromosome, thus even enabling using reads that cover only individual variants. We demonstrate HapTree-X’s feasibility on in-house sequenced Genome in a Bottle RNA-seq and various whole exome, genome, and 10X Genomics datasets. HapTree-X produces more complete phases (up to 25%), even in clinically important genes, and phases more variants than other methods while maintaining similar or higher accuracy and being up to 10×  faster than other tools. The advantage of HapTree-X’s ability to use multiple lines of evidence, as well as to phase polyploid genomes in a single integrative framework, substantially grows as the amount of diverse data increases.

ASEP: Gene-based detection of allele-specific expression across individuals in a population by RNA sequencing.

Allele-specific expression (ASE) analysis, which quantifies the relative expression of two alleles in a diploid individual, is a powerful tool for identifying cis-regulated gene expression variations that underlie phenotypic differences among individuals. Existing methods for gene-level ASE detection analyze one individual at a time, therefore failing to account for shared information across individuals. Failure to accommodate such shared information not only reduces power, but also makes it difficult to interpret results across individuals. However, when only RNA sequencing (RNA-seq) data are available, ASE detection across individuals is challenging because the data often include individuals that are either heterozygous or homozygous for the unobserved cis-regulatory SNP, leading to sample heterogeneity as only those heterozygous individuals are informative for ASE, whereas those homozygous individuals have balanced expression. To simultaneously model multi-individual information and account for such heterogeneity, we developed ASEP, a mixture model with subject-specific random effect to account for multi-SNP correlations within the same gene. ASEP only requires RNA-seq data, and is able to detect gene-level ASE under one condition and differential ASE between two conditions (e.g., pre- versus post-treatment). Extensive simulations demonstrated the convincing performance of ASEP under a wide range of scenarios. We applied ASEP to a human kidney RNA-seq dataset, identified ASE genes and validated our results with two published eQTL studies. We further applied ASEP to a human macrophage RNA-seq dataset, identified genes showing evidence of differential ASE between M0 and M1 macrophages, and confirmed our findings by results from cardiometabolic trait-relevant genome-wide association studies. To the best of our knowledge, ASEP is the first method for gene-level ASE detection at the population level that only requires the use of RNA-seq data. With the growing adoption of RNA-seq, we believe ASEP will be well-suited for various ASE studies for human diseases.

Abstract 98: Harnessing Novel Gene Expression Analyses To Identify Drivers Of Ear Regeneration In MRL Mice

Purpose: The “super-healing” MRL/MpJ mouse is uniquely able to regenerate ear punch wounds, while dorsal wounds heal in a typical scarring fashion. Identification of the genes responsible for ear regeneration in these mice could provide insights toward novel anti-scarring therapies. However, previous methods for identifying the genetic basis of species-/strain-specific traits are imprecise and extremely resource-intensive. Recent studies have revealed that, by outbreeding a strain with a trait of interest (here, regenerative ear healing), assessing the hybrid offspring for spatial/tissue-specific differences in relative gene expression from the two parent alleles (“differential allele-specific expression”) can reveal genes underlying that phenotype. We hypothesized that differential allele-specific expression could represent a novel approach to identify genes driving regenerative ear healing in MRL mice. Methods: MRL/MpJ (MRL) and CAST/EiJ (CAST) mouse genome sequences were obtained. Variant calling was performed to identify sites where the sequences differed. MRL and CAST mice were cross-bred, F1 offspring underwent dorsal excisional and ear punch wounding (standard protocols), and wounds were harvested on POD7. Following incubation in ammonium thiocyanate, epidermis and underlying cartilage (for ear wounds) were removed under a surgical microscope. Dorsal and ear wound dermis were digested, and inflammatory/immune cells (CD45+; chosen because prior studies have implicated altered wound inflammatory response in MRL regeneration) were isolated by FACS and subjected to bulk RNA-sequencing. Reads were mapped to the MRL or CAST genome based on strain-specific variants, and genes with significantly differing expression (Benjamini-Hochberg method, FDR=5%) from the MRL vs. CAST allele and from dorsal vs. ear wounds were identified. Results: Ear wounds in MRLxCAST hybrid mice exhibited grossly apparent regeneration over time, consistent with prior reports. Roughly equal proportions of RNA-seq reads mapped to MRL vs. CAST alleles, indicating minimal mapping bias. Principal component analysis demonstrated that reads clustered by parent allele (MRL vs. CAST) and wound site (ear vs. dorsum), with good reproducibility between biological replicates (n=4). Importantly, the number of genes with significantly different expression from MRL vs. CAST allele was substantially (>10-fold) greater in the ear than the dorsum (2531 vs. 159 unique genes. The majority of allele-specific expression in the dorsum was also observed in the ear, overlapping region; 878 shared genes vs. 159 unique to dorsum). Conclusions: In MRL hybrid mice, inflammatory cells exhibit markedly greater allele-specific gene expression in ear wounds (regenerative) than dorsal wounds (scarring). This finding would be consistent with MRL-specific gene cis-regulatory activity leading to the improved healing phenotype seen in the ear. Our results suggest that probing differential allele-specific gene expression may represent a powerful computational approach for identifying the genetic basis of complex traits in mammals. We will examine inflammatory and other cell types (e.g., fibroblasts) in order to identify specific genes with differential allele-specific expression in regenerative versus scarring wounds, with the ultimate goal of identifying novel targets for preventing scarring.

Differential Allele-Specific Expression Uncovers Breast Cancer Genes Dysregulated by Cis Noncoding Mutations.

Identifying cancer-relevant mutations in noncoding regions is challenging due to the large numbers of such mutations, their low levels of recurrence, and difficulties in interpreting their functional impact. To uncover genes that are dysregulated due to somatic mutations in cis, we build upon the concept of differential allele-specific expression (ASE) and introduce methods to identify genes within an individual's cancer whose ASE differs from what is found in matched normal tissue. When applied to breast cancer tumor samples, our methods detect the known allele-specific effects of copy number variation and nonsense-mediated decay. Further, genes that are found to recurrently exhibit differential ASE across samples are cancer relevant. Genes with cis mutations are enriched for differential ASE, and we find 147 potentially functional noncoding mutations cis to genes that exhibit significant differential ASE. We conclude that differential ASE is a promising means for discovering gene dysregulation due to cis noncoding mutations.

Abstract 1584: Transcriptome analysis links immune genes allelic expression imbalances to lung cancer

Abstract Background: Genome-wide association study (GWAS) have identified over 45 susceptibility loci for lung cancer; many studies including our own group, have focused on low-frequency and rare coding variants using fine mapping and exome sequencing. This strategy, however, has met with limited success as about 90% of GWAS hits are noncoding and act primarily through altering transcriptional regulation in an allele-specific manner. The RNA-Seq based allele-specific expression (ASE) analysis affords an innovative approach to study preferential expression of an allele in direct relationship to its genotype, providing information on cis-regulatory effects for the expression of putative genes. However currently, there are no lung cancer studies that have rigorously evaluated the ASE variation in lung tumor and adjacent tissues. Methods: Leveraging The Cancer Genome Atlas (TCGA) resource, we performed transcriptomic-wide ASE analysis using existing RNA-Seq datasets of paired tumor and adjacent tissues from 54 lung adenocarcinoma patients. We first quantified the RNA read counts of Referent and Alternate alleles of heterozygous variants, then evaluated the allelic imbalance on a per-sample basis using Beta-binomial test, and explored the differential ASE between tumor and adjacent tissues using paired Wilcoxon test. Functional regulatory consequences were generated from Ensembl Variant Effect Predictor. Results: We identified total 208 significant ASEs, including 35 tissue-specific (only in tumor or only in adjacent), 28 sharing, and 145 differential variants. Of the 208 candidates, 41 were from the human leukocyte antigen (HLA) locus (primary DQA2, DQB1, DRB1, H and J), 26 were from the immunoglobulin (IG) superfamily (primary IGH, IGL, IGK and F11R). About 80% candidates were noncoding (mostly in 5’ and 3’ untranslated regions) and with regulatory features (21 promoter, seven enhancer, seven open chromatin region, two induce nonsense-mediated mRNA decay, one CTCF-binding site, and one transcription factor binding site). Other top genes included MDM2, APOL1, and CTSB. Pathway analyses revealed 27 genes involved in immune response pathway, and 12 genes involved in HLA antigen processing and presentation pathway. Conclusion: This study is the first transcriptomics ASE analysis in lung adenocarcinoma. The key somatic cis-regulatory ASE variants identified from this study, especially immunogenic allelic variations from HLA and IG genes, could be used for identifying high-risk individuals for targeted lung cancer checkpoint blockade and related immunotherapies. Citation Format: Yanhong Liu, Spiridon Tsavachidis, Farrah Kheradmand, Margaret R. Spitz, Chris Amos. Transcriptome analysis links immune genes allelic expression imbalances to lung cancer [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2019; 2019 Mar 29-Apr 3; Atlanta, GA. Philadelphia (PA): AACR; Cancer Res 2019;79(13 Suppl):Abstract nr 1584.

Fully Bayesian analysis of allele-specific RNA-seq data.

Diploid organisms have two copies of each gene, called alleles, that can be separately transcribed. The RNA abundance associated to any particular allele is known as allele-specific expression (ASE). When two alleles have polymorphisms in transcribed regions, ASE can be studied using RNA-seq read count data. ASE has characteristics different from the regular RNA-seq expression: ASE cannot be assessed for every gene, measures of ASE can be biased towards one of the alleles (reference allele), and ASE provides two measures of expression for a single gene for each biological samples with leads to additional complications for single-gene models. We present statistical methods for modeling ASE and detecting genes with differential allelic expression. We propose a hierarchical, overdispersed, count regression model to deal with ASE counts. The model accommodates gene-specific overdispersion, has an internal measure of the reference allele bias, and uses random effects to model the gene-specific regression parameters. Fully Bayesian inference is obtained using the fbseq package that implements a parallel strategy to make the computational times reasonable. Simulation and real data analysis suggest the proposed model is a practical and powerful tool for the study of differential ASE.

Hypoxic Response-Dependent Genetic Regulation Revealed By Allele-Specific Expression in Reticulocytes of Chuvash Polycythemia

Homozygosity for the VHLR200W mutation in Chuvash polycythemia (CP) leads to decreased degradation of the α subunits of hypoxia inducible factor (HIF)-1 and HIF-2 by the hypomorphic variant of VHL, the principal negative regulator of HIFs. The constitutively activated HIFs directly regulate the transcription of a suite of hypoxic responsible genes, including the principal regulators of erythropoiesis, vessel development, and glycolytic metabolism, which further trigger a downstream cascade of gene expression. Besides these transcriptional factors, cis acting elements play an important role in the hypoxic gene regulatory network. To assess the extent of cis regulatory variation in hypoxic gene expression, we compared allele-specific expression (ASE) in purified reticulocytes between VHLR200W homozygote individuals and age- and gender-matched wild type control individuals living at the same altitude of ~200 meters from the Chuvash population.Cell fractions of reticulocytes were purified from 17 VHLR200W homozygotes and 13 wild type individuals. Total RNA was extracted, depleted of ribosomal RNA and hemoglobin transcripts, and reverse transcribed. Strand-specific libraries were constructed for 125 bp paired-end sequencing to 30-45 million read pairs per sample using Illumina HiSeq 2500 platform. The samples were collected and processed in three batches across two years, with VHL genotype randomized in each batch. The sequencing data were mapped to human reference genome and analyzed for differential expression and differential ASE between VHLR200W homozygotes and wild type individuals.At 5% false discovery rate (FDR, i.e., <5 false positives in 100 detected genes), 1,267 genes were differentially expressed with more than 1.2-fold change in CP patients, 703 elevated and 564 decreased. Genes up-regulated in CP were enriched (fold enrichment >5, FDR <0.05) in REACTOME pathways of epigenetic remodeling (Packaging of telomere ends, DNA methylation, HDACs deacetylate histones, PRC2 methylates histones and DNA, Deposition of new CENPA-containing nucleosomes at the centromere, HATs acetylate histones) and oxidative stress induced senescence (DNA damage/telomere stress induced senescence, Senescence-associated secretory phenotype, Oxidative stress induced senescence). Genes decreased in CP were enriched in REACTOME pathways of cell cycle (E2F-enabled inhibition of pre-replication complex formation, Nuclear pore complex disassembly, SUMOylation of DNA replication proteins) and DNA damage repair (Activation of ATR in response to replication stress, SUMOylation of DNA damage response and repair proteins).ASE was analyzed between CP and wild type individuals to assess hypoxic response-dependent genetic effects on gene expression. For the 1,267 genes differentially expressed in the CP, we selected genes containing exonic SNPs with heterozygous alleles for ASE analysis. With a null hypothesis of no cis acting regulation on the gene expression, both alleles are expected to be expressed at the same level, whereas allelic imbalance indicates linked cis regulation. At a given bi-allelic SNP, individuals with ≥2 read counts covering each of the reference and alternative alleles and with ≥20 total counts were included in the analysis. Exonic SNPs with at least one individual in each of the CP and wild type group were further selected to test for differential ASE between the CP and wild type groups, using a generalized linear model. A total of 147 genes passed the filtering and were analyzed, among which 32 were detected to have significant CP-dependent ASE at 5% FDR. Some of these genes may have important roles in hypoxic responses in CP reticulocytes, for example NEIL3, encoding a DNA glycosylase that initiates the first step in base excision repair by cleaving bases damaged by reactive oxygen species, and STOM, encoding an integral membrane protein that localizes to the cell membrane of red blood cells, loss of which is associated with hereditary stomatocytosis.Our study reveals plethora of gene expression changes in CP reticulocytes compared to wild type controls, among which 22% could be regulated by hypoxic response-specific cis genetic variations. These observations indicate the prominence of cis elements in hypoxic response, for which substantial inter-individual differences exist even among a relatively isolated population. DisclosuresGordeuk:Emmaus Life Sciences: Consultancy.

Abstract 1884: The breast cancer gene CCDC170 regulates the Golgi-associated microtubule network and directional cell migration

Abstract Background: Breast cancer (BCa) is a genetically heterogeneous disease and many genes contributing to BCa risk remain to be identified. Genome-Wide Association Studies (GWAS) and subsequent fine-mapping studies (&amp;gt;50) have strongly implicated genetic alterations at the CCDC170/C6ORF97-ESR1 locus (6q25.1) as being associated with the risk of BCa. ESR1, encoding the estrogen receptor, might be a more obvious candidate for influencing risk. Surprisingly, our analysis using genome-wide differential allele-specific expression (DASE), an indicator for BCa susceptibility, suggested that the genetic alterations of CCDC170, but not ESR1, accounts for GWAS-associated BCa risk at this locus. CCDC170 is a coiled-coiled domain-containing protein of poorly understood function. BCa-specific truncation and missense mutations in CCDC170 also have been detected, with the truncations being implicated in driving Luminal B subtype BCa. Together these findings demonstrate that the CCDC170 gene is involved in BCa, but the underlying molecular mechanisms for its role in tumorigenesis are unknown. Experimental designs and results: By using the approaches of confocal microscopy and cell imaging analysis, here we report for the first time that CCDC170 is associated with the Golgi apparatus and perinuclear microtubules (MTs), and support a role for CCDC170 in the Golgi-associated microtubule network. We have shown that overexpression of CCDC170 triggers Golgi reorganization and stabilizes Golgi-associated MTs, accompanied by dramatically increased acetylation of α-tubulin that is driven by the acetyltransferase ATAT1. The Golgi-associated MT network has been proposed to regulate cell polarity and migration. In support of this concept, we have shown that CRISPR knockout increases, and overexpression of CCDC170 decreases, BCa directional cell migration in vitro. We also found that the BCa-specific truncations result in mislocalization of CCDC170 and/or diminished stability of Golgi-associated MTs. Lastly, we identified candidate CCDC170 functional binding partners (e.g. MAP4) that are consistent with its localization and proposed function. These partners may serve to mediate the acetylation and stabilization of MTs. Conclusions: Taken together, our findings demonstrate that CCDC170 plays an essential role in Golgi-associated MT organization and stabilization, and provide a mechanism for how perturbations in CCDC170 could alter Golgi-mediated cell polarity, and thereby drive BCa and other abnormalities. This work was partially supported by the Susan G. Komen for the Cure (KG100274), NCI (CA186853), and Eileen Stein Jacoby Fund. Citation Format: Pengtao Jiang, Yueran Li, Andrey Poleshko, Valentina Medvedeva, Natalia Baulina, Yongchao Zhang, Yan Zhou, Carolyn M. Slater, Trinity Pellegrin, Jason Wasserman, Michael Lindy, Mary Daly, Richard A. Katz, Xiaowei Chen. The breast cancer gene CCDC170 regulates the Golgi-associated microtubule network and directional cell migration [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2017; 2017 Apr 1-5; Washington, DC. Philadelphia (PA): AACR; Cancer Res 2017;77(13 Suppl):Abstract nr 1884. doi:10.1158/1538-7445.AM2017-1884

The Protein Encoded by the CCDC170 Breast Cancer Gene Functions to Organize the Golgi-Microtubule Network

Genome-Wide Association Studies (GWAS) and subsequent fine-mapping studies (>50) have implicated single nucleotide polymorphisms (SNPs) located at the CCDC170/C6ORF97-ESR1 locus (6q25.1) as being associated with the risk of breast cancer. Surprisingly, our analysis using genome-wide differential allele-specific expression (DASE), an indicator for breast cancer susceptibility, suggested that the genetic alterations of CCDC170, but not ESR1, account for GWAS-associated breast cancer risk at this locus. Breast cancer-associated CCDC170 nonsense mutations and rearrangements have also been detected, with the latter being specifically implicated in driving breast cancer. Here we report that the wild type CCDC170 protein localizes to the region of the Golgi apparatus and binds Golgi-associated microtubules (MTs), and that breast cancer-linked truncations of CCDC170 result in loss of Golgi localization. Overexpression of wild type CCDC170 triggers Golgi reorganization, and enhances Golgi-associated MT stabilization and acetyltransferase ATAT1-dependent α-tubulin acetylation. Golgi-derived MTs regulate cellular polarity and motility, and we provide evidence that dysregulation of CCDC170 affects polarized cell migration. Taken together, our findings demonstrate that CCDC170 plays an essential role in Golgi-associated MT organization and stabilization, and implicate a mechanism for how perturbations in the CCDC170 gene may contribute to the hallmark changes in cell polarity and motility seen in breast cancer.

Allele-specific expression in the human heart and its application to postoperative atrial fibrillation and myocardial ischemia.

BackgroundAllele-specific expression (ASE) is differential expression of each of the two chromosomal alleles of an autosomal gene. We assessed ASE patterns in the human left atrium (LA, n = 62) and paired samples from the left ventricle (LV, n = 76) before and after ischemia, and tested the utility of differential ASE to identify genes associated with postoperative atrial fibrillation (poAF) and myocardial ischemia.MethodsFollowing genotyping from whole blood and whole-genome sequencing of LA and LV samples, we called ASE using sequences overlapping heterozygous SNPs using rigorous quality control to minimize false ASE calling. ASE patterns were compared between cardiac chambers and with a validation cohort from cadaveric tissue. ASE patterns in the LA were compared between patients who had poAF and those who did not. Changes in ASE in the LV were compared between paired baseline and post-ischemia samples.ResultsASE was found for 3404 (5.1%) and 8642 (4.0%) of SNPs analyzed in the LA and LV, respectively. Out of 6157 SNPs with ASE analyzed in both chambers, 2078 had evidence of ASE in both LA and LV (p < 0.0001). The SNP with the greatest ASE difference in the LA of patients with and without postoperative atrial fibrillation was within the gelsolin (GSN) gene, previously associated with atrial fibrillation in mice. The genes with differential ASE in poAF were enriched for myocardial structure genes, indicating the importance of atrial remodeling in the pathophysiology of AF. The greatest change in ASE between paired post-ischemic and baseline samples of the LV was in the zinc finger and homeodomain protein 2 (ZHX2) gene, a modulator of plasma lipids. Genes with differential ASE in ischemia were enriched in the ubiquitin ligase complex pathway associated with the ischemia-reperfusion response.ConclusionsOur results establish a pattern of ASE in the human heart, with a high degree of shared ASE between cardiac chambers as well as chamber-specific ASE. Furthermore, ASE analysis can be used to identify novel genes associated with (poAF) and myocardial ischemia.Electronic supplementary materialThe online version of this article (doi:10.1186/s13073-016-0381-1) contains supplementary material, which is available to authorized users.

Allelic expression imbalance polymorphisms in susceptibility chromosome regions and the risk and survival of breast cancer.

Allelic expression imbalance (AEI) has been applied to indicate potential function of genetic variants. Combining earlier results from global differential allele-specific expression analysis and genome wide association studies (GWASs), we select the single nuclear polymorphisms (SNPs) exhibiting AEI phenomenon located in breast cancer susceptibility chromosome regions, and evaluate their associations with breast cancer risk and survival. We examined the genotypes of 10 AEI SNPs in 1551 incident breast cancer cases and 1605 age-frequency matched controls from Guangzhou, China. In total, 1168 cases were followed up. MUC16 rs2591592 (AT/AA vs. TT) was associated with an increased risk of premenopausal breast cancer (OR [95%CI]: 1.30 [1.07, 1.57]); SLAMF1 rs1061217 (CT/TT vs. CC) decreased the risk of breast cancer among overweight women (OR [95%CI]: 0.74 [0.57, 0.96]) but increased the risk among normal-weight women (OR [95%CI]: 1.15 [1.01, 1.39]); ZNF331 rs8109631 (AG/AA vs. GG) and CHRAC1 rs10216653 (GC/GG vs. CC) were associated with progression free survival among breast cancer patients with negative ER/PR status and higher clinical stage (HRs [95%CIs]: 2.39 [1.14, 5.00], 1.85 [1.03, 3.32], and 0.49 [0.30, 0.80], respectively). ZNF331 rs8109631 and CHRAC1 rs10216653 were further found to represent several functional SNPs through bioinformatic analysis. In conclusion, our findings demonstrated suggestive associations of AEI polymorphisms with breast cancer risk (MUC16 rs2591592 and SLAMF1 rs1061217) and prognosis (ZNF331 rs8109631 and CHRAC1 rs10216653). © 2016 Wiley Periodicals, Inc.

A GWAS SNP for Schizophrenia Is Linked to the Internal MIR137 Promoter and Supports Differential Allele-Specific Expression

Single nucleotide polymorphisms (SNPs) within the MIR137 gene locus have been shown to confer risk for schizophrenia through genome-wide association studies (GWAS). The expression levels of microRNA-137 (miR-137) and its validated gene targets have also been shown to be disrupted in several neuropsychiatric conditions, including schizophrenia. Regulation of miR-137 expression is thus imperative for normal neuronal functioning. We previously characterized an internal promoter domain within the MIR137 gene that contained a variable number tandem repeat (VNTR) polymorphism and could alter the in vitro levels of miR-137 in a stimulus-induced and allele-specific manner. We now demonstrate that haplotype tagging-SNP analysis linked the rs1625579 GWAS SNP for schizophrenia to this internal MIR137 promoter through a proxy SNP rs2660304 located at this domain. We postulated that the rs2660304 promoter SNP may act as predisposing factor for schizophrenia through altering the levels of miR-137 expression in a genotype-dependent manner. Reporter gene analysis of the internal MIR137 promoter containing the common VNTR variant demonstrated genotype-dependent differences in promoter activity with respect to rs2660304. In line with previous reports, the major allele of the rs2660304 proxy SNP, which has previously been linked with schizophrenia risk through genetic association, resulted in downregulation of reporter gene expression in a tissue culture model. The genetic influence of the rs2660304 proxy SNP on the transcriptional activity of the internal MIR137 promoter, and thus the levels of miR-137 expression, therefore offers a distinct regulatory mechanism to explain the functional significance of the rs1625579 GWAS SNP for schizophrenia risk.

Identifying breast cancer risk loci by global differential allele-specific expression (DASE) analysis in mammary epithelial transcriptome

BackgroundThe significant mortality associated with breast cancer (BCa) suggests a need to improve current research strategies to identify new genes that predispose women to breast cancer. Differential allele-specific expression (DASE) has been shown to contribute to phenotypic variables in humans and recently to the pathogenesis of cancer. We previously reported that nonsense-mediated mRNA decay (NMD) could lead to DASE of BRCA1/2, which is associated with elevated susceptibility to breast cancer. In addition to truncation mutations, multiple genetic and epigenetic factors can contribute to DASE, and we propose that DASE is a functional index for cis-acting regulatory variants and pathogenic mutations, and that global analysis of DASE in breast cancer precursor tissues can be used to identify novel causative alleles for breast cancer susceptibility.ResultsTo test our hypothesis, we employed the Illumina® Omni1-Quad BeadChip in paired genomic DNA (gDNA) and double-stranded cDNA (ds-cDNA) samples prepared from eight BCa patient-derived normal mammary epithelial lines (HMEC). We filtered original array data according to heterozygous genotype calls and calculated DASE values using the Log ratio of cDNA allele intensity, which was normalized to the corresponding gDNA. We developed two statistical methods, SNP- and gene-based approaches, which allowed us to identify a list of 60 candidate DASE loci (DASE ≥ 2.00, P ≤ 0.01, FDR ≤ 0.05) by both methods. Ingenuity Pathway Analysis of DASE loci revealed one major breast cancer-relevant interaction network, which includes two known cancer causative genes, ZNF331 (DASE = 2.31, P = 0.0018, FDR = 0.040) and USP6 (DASE = 4.80, P = 0.0013, FDR = 0.013), and a breast cancer causative gene, DMBT1 (DASE=2.03, P = 0.0017, FDR = 0.014). Sequence analysis of a 5′ RACE product of DMBT1 demonstrated that rs2981745, a putative breast cancer risk locus, appears to be one of the causal variants leading to DASE in DMBT1.ConclusionsOur study demonstrated for the first time that global DASE analysis is a powerful new approach to identify breast cancer risk allele(s).

Abstract LB-407: Global differential allele-specific expression (DASE) analysis in mammary epithelial transcriptome

Abstract The significant mortality associated with breast cancer (BCa) suggests a clear need to improve current research strategies to identify new genes that predispose women to breast cancer. Differential allele-specific expression (DASE) has been shown to contribute to phenotypic variables in human and recently to the pathogenesis of cancer. We previously reported that nonsense-mediated mRNA decay (NMD) could lead to DASE of BRCA1/2 which is associated with elevated susceptibility to BCa. In addition to truncation mutations, multiple genetic and epigenetic factors, such as mutations in promoter regions, alteration in microRNA binding sites on 3′UTRs and changes in chromatin methylation status can contribute to DASE. Therefore, we hypothesize that DASE is a functional index for cis-acting variants and pathogenic mutations and global profiling of DASE in BCa precursor tissues can be used to identify novel causative alleles for BCa susceptibility. To test our hypothesis, we employed the Illumina® Omni1-Quad BeadChip in paired genomic DNA (gDNA) and double-stranded cDNA (ds-cDNA) samples prepared from of eight BCa patient-derived normal mammary epithelial lines (HMEC). We filtered original array data according to heterozygous genotype calls from 86,267 female transcriptome SNPs covered by the BeadChip and calculated DASE values using the Log ratio of cDNA allele intensity, which was normalized to the corresponding gDNA intensity to limit the potential technical variation cause by the design of oligo probes. We have developed two statistical methods, SNP- and gene-based approaches, which allowed us to identify 90 and 143 genes with a significant DASE&amp;gt;2 (P &amp;lt;= 0.01, FDR &amp;lt;= 0.05), respectively. Importantly, 60 candidate loci are identified by both approaches. Ingenuity Pathway Analysis of the 34 coding genes in top 60 candidates revealed two major interaction networks, which involves a variety of biological functions, including cell death and cell-to-cell signaling. One of the networks includes known cancer causative gene ZNF331 (DASE = 2.31, P = 0.0018, FDR = 0.040), USP6 (DASE = 4.80, P = 0.0013, FDR = 0.013) and breast cancer causative gene DMBT1 (DASE=2.03, P = 0.0017, FDR = 0.014). Sequencing of the 5′UTR of DMBT1 identified a genetic variant, rs2981745, which has been reported previously as a causative variant for breast cancer risk in DMBT1 locus. Further sequencing of 5′ RACE product of DMBT1 confirmed rs2981745 as the causal variant leading to DASE in DMBT1. This finding indicates that abnormalities in phenotypes of DASE reflect pathogenic genotypes. In summary, our study clearly demonstrated that global DASE profiling is a powerful new approach to identify BCa risk alleles. This work was kindly supported by the Susan G. Komen for the Cure (KG100274) and Eileen Stein Jacoby Fund. Citation Format: {Authors}. {Abstract title} [abstract]. In: Proceedings of the 103rd Annual Meeting of the American Association for Cancer Research; 2012 Mar 31-Apr 4; Chicago, IL. Philadelphia (PA): AACR; Cancer Res 2012;72(8 Suppl):Abstract nr LB-407. doi:1538-7445.AM2012-LB-407

Quantitative measurement of allele-specific protein expression in a diploid yeast hybrid by LC-MS.

Understanding the genetic basis of gene regulatory variation is a key goal of evolutionary and medical genetics. Regulatory variation can act in an allele-specific manner (cis-acting) or it can affect both alleles of a gene (trans-acting). Differential allele-specific expression (ASE), in which the expression of one allele differs from another in a diploid, implies the presence of cis-acting regulatory variation. While microarrays and high-throughput sequencing have enabled genome-wide measurements of transcriptional ASE, methods for measurement of protein ASE (pASE) have lagged far behind. We describe a flexible, accurate, and scalable strategy for measurement of pASE by liquid chromatography-coupled mass spectrometry (LC-MS). We apply this approach to a hybrid between the yeast species Saccharomyces cerevisiae and Saccharomyces bayanus. Our results provide the first analysis of the relative contribution of cis-acting and trans-acting regulatory differences to protein expression divergence between yeast species.

Differential Allelic Expression of Dopamine D1 Receptor Gene ( DRD1) Is Modulated by microRNA miR-504

Previously, we reported that dopamine D1 receptor gene (DRD1) is associated with nicotine dependence (ND) and demonstrated that two alleles (A and G) of polymorphism rs686 in the 3'-untranslated region (3'UTR) of DRD1 are expressed differentially. However, the mechanism underlying the differential expression remains to be determined. We hypothesize that it is caused by miRNA targeting. We first used the MicroInspector algorithm to identify microRNAs (miRNAs) potentially targeting the rs686 polymorphism in the DRD1 3'UTR and then employed a luciferase reporter assay combined with site-directed mutagenesis to test the predicted miRNA targeting. We also examined the miRNA targeting of DRD1 with a gene expression assay. Of two miRNAs predicted by computational analyses, we found that miR-504, not miR-296, upregulated reporter luciferase activity and increased DRD1 expression by targeting the DRD1 3' UTR, whereas inhibition of miR-504, not miR-296, had the opposite effect. Furthermore, we showed that the direct binding of miR-504 to the DRD1 3'UTR, verified by site-directed mutagenesis, causes a significant expression difference between the two alleles. miR-504 up-regulates DRD1 expression by direct binding to the 3'UTR, which leads to differential allele-specific expression of DRD1.

Differential allele expression of host defense genes, pulmonary surfactant protein-A and Osteopontin, in rat

Differential allele-specific expression has been observed in several genes involved in immunity. SP-A and OPN play a role in innate host defense. To determine whether SP-A and OPN are subject to differential allele-specific regulation, we investigated their gene or allele-specific expression in various tissues. The results showed: (1) Tissue-specific expression with high levels in lung ( SP-A) and kidney ( OPN). (2) Differences in allele-specific expression among individuals and tissues. SP-A showed an exclusively balanced biallelic expression (BB) in lung, but both BB and imbalanced biallelic (IB) expression in colon. Allele expression of OPN was more heterogeneous, e.g. in colon BB (22%), IB (64%), and monoallelic expression (MO) (14%). (3) Differential allele-specific expression was observed in all tissues studied ( OPN) or in all extrapulmonary tissues ( SP-A). (4) Family studies indicated that inheritable factor(s) may be involved in the regulation of allele-specific expression. (5) Analysis of co-expression of gene-specific alleles from double heterozygous rats revealed lack of coordinate allele expression among SP-A, SP-D, and OPN. We conclude that allele-specific expression occurs among genes of innate host defense. This may yet provide another level of regulatory complexity for molecules involved in the first line of defense.